----------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_table1.log
  log type:  text
 opened on:  23 Jun 2022, 23:05:39

. 
.         use "$Rep_smokelabor/1_build/regdata/county_quarter.dta", clear 

.         
.         local tbl_settings_log format(%6.3f) parentheses(stderr) asterisk()

.         local tbl_settings_pmil format(%6.1f) parentheses(stderr) asterisk()

.         
.         ** panel a 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_dens [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      89.59
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7883
                                                  Adj R-squared   =     0.7655
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0454
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.3422

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_dens |    .053042   .0056039     9.47   0.000     .0420507    .0640333
       _cons |   9.109334   .0449067   202.85   0.000     9.021256    9.197413
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", table(pm25, `tbl_settings_log') 
> replace
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_a.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_dens [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      27.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2426
                                                  Adj R-squared   =     0.1754
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0032
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   409.7888

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_dens |  -4.102343   .7789493    -5.27   0.000    -5.629789   -2.574896
       _cons |   168.1804   5.860071    28.70   0.000     156.6894    179.6715
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", table(earnings, `tbl_settings_lo
> g') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_a.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_dens [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stq
> tros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      14.00
Statistics robust to heteroskedasticity           Prob > F        =     0.0002
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4876
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0005
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 18019.4830

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_dens |  -74.98713   20.04266    -3.74   0.000    -114.2889   -35.68535
       _cons |  -1283.021   150.5435    -8.52   0.000    -1578.223   -987.8193
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", table(employment, `tbl_settings_
> pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_a.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_dens [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   2547) =       6.87
Statistics robust to heteroskedasticity           Prob > F        =     0.0088
                                                  R-squared       =     0.2402
                                                  Adj R-squared   =     0.1734
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0002
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11152.3131

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_dens |  -25.28231   9.643659    -2.62   0.009    -44.19253   -6.372103
       _cons |  -1963.761   80.53286   -24.38   0.000    -2121.677   -1805.844
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", table(lfp, `tbl_settings_pmil') 
> append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_a.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_dens_coef   0.053***   -4.102***   -75.0***   -25.3*** |
  2. | hms_dens_stderr    (0.006)     (0.779)     (20.0)      (9.6) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_a.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_a.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel b, first row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 g*_tavg g*_wdir hms_deep [aw=seer_pop], a(g_ppt g_wspd fe_countyqtroy fe_styr) vce(cluster countyfi
> p fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 16 iterations)

HDFE Linear regression                            Number of obs   =     75,181
Absorbing 4 HDFE groups                           F(  16,   1685) =      23.50
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.8117
                                                  Adj R-squared   =     0.7913
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0756
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.2660

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     g1_tavg |    .023198   .0102952     2.25   0.024     .0030054    .0433907
     g2_tavg |   .0795833   .0106675     7.46   0.000     .0586604    .1005062
     g3_tavg |   .0106669   .0079935     1.33   0.182    -.0050112    .0263451
     g4_tavg |   .0110791   .0068176     1.63   0.104    -.0022928    .0244509
     g5_tavg |  -.0132983   .0089252    -1.49   0.136    -.0308038    .0042073
     g6_tavg |  -.0166333   .0049096    -3.39   0.001    -.0262628   -.0070037
     g7_tavg |    .006354   .0051689     1.23   0.219    -.0037842    .0164922
     g8_tavg |   .0151775   .0036699     4.14   0.000     .0079795    .0223755
     g9_tavg |    .010281   .0054924     1.87   0.061    -.0004917    .0210537
     g1_wdir |    .093787    .035986     2.61   0.009     .0232049     .164369
     g2_wdir |   .0994946   .0365208     2.72   0.007     .0278637    .1711254
     g3_wdir |   .0980377    .036355     2.70   0.007      .026732    .1693435
     g4_wdir |   .1145532   .0378714     3.02   0.003     .0402734    .1888331
     g5_wdir |   .1321042   .0366562     3.60   0.000     .0602076    .2040008
     g6_wdir |   .1073913   .0382329     2.81   0.005     .0324024    .1823802
    hms_deep |    .043002   .0063955     6.72   0.000     .0304581    .0555459
       _cons |  -.6808793   3.301562    -0.21   0.837    -7.156474    5.794716
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
          g_ppt |        10           0          10     |
         g_wspd |        10           1           9     |
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           1         636    ?|
--------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", table(pm25, `tbl_setti
> ngs_log') replace
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll g*_tavg g*_wdir hms_deep [aw=seer_pop], a(g_ppt g_wspd fe_countyqtroy fe_styr) vce(clus
> ter countyfip fe_stqtros)
(MWFE estimator converged in 15 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =    146,268
Absorbing 4 HDFE groups                           F(  16,   2519) =     451.38
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2502
                                                  Adj R-squared   =     0.1831
Number of clusters (countyfip) =      2,837       Within R-sq.    =     0.0088
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   407.6070

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     g1_tavg |  -6.653556   1.527288    -4.36   0.000    -9.648424   -3.658688
     g2_tavg |  -4.902943   1.901721    -2.58   0.010    -8.632039   -1.173848
     g3_tavg |  -4.350799   1.039289    -4.19   0.000    -6.388748    -2.31285
     g4_tavg |  -4.413547   .8908674    -4.95   0.000    -6.160454   -2.666639
     g5_tavg |  -1.037041   .7113886    -1.46   0.145    -2.432007     .357925
     g6_tavg |  -1.161983   .8117052    -1.43   0.152     -2.75366    .4296952
     g7_tavg |  -.2447154   .6239046    -0.39   0.695    -1.468134    .9787029
     g8_tavg |    .786563   .6272744     1.25   0.210    -.4434632    2.016589
     g9_tavg |   4.350432   .8144569     5.34   0.000     2.753359    5.947506
     g1_wdir |  -10.80851   2.811962    -3.84   0.000     -16.3225   -5.294516
     g2_wdir |  -9.998193   2.830347    -3.53   0.000    -15.54824   -4.448147
     g3_wdir |  -12.16843   2.731599    -4.45   0.000    -17.52484   -6.812021
     g4_wdir |  -11.40311   2.734136    -4.17   0.000     -16.7645    -6.04173
     g5_wdir |  -12.02761    2.99392    -4.02   0.000    -17.89841   -6.156818
     g6_wdir |  -13.01825   3.018388    -4.31   0.000    -18.93702   -7.099469
    hms_deep |  -4.950374   .7911605    -6.26   0.000    -6.501766   -3.398983
       _cons |   1255.918   248.7166     5.05   0.000     768.2082    1743.628
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
          g_ppt |        10           0          10     |
         g_wspd |        10           1           9     |
 fe_countyqtroy |     11348       11348           0    *|
        fe_styr |       631           1         630    ?|
--------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", table(earnings, `tbl_s
> ettings_log') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal g*_tavg g*_wdir hms_deep [aw=seer_pop16plus], a(g_ppt g_wspd fe_countyqtroy fe_styr)
>  vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 15 iterations)

HDFE Linear regression                            Number of obs   =    146,268
Absorbing 4 HDFE groups                           F(  16,   2519) =      14.92
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.5415
                                                  Adj R-squared   =     0.5005
Number of clusters (countyfip) =      2,837       Within R-sq.    =     0.0090
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 17697.3786

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     g1_tavg |  -206.9054   72.96449    -2.84   0.005    -349.9819   -63.82886
     g2_tavg |   -227.943   67.19988    -3.39   0.001    -359.7156    -96.1703
     g3_tavg |  -229.1851   53.05183    -4.32   0.000    -333.2148   -125.1555
     g4_tavg |   -218.585   54.20026    -4.03   0.000    -324.8666   -112.3034
     g5_tavg |  -74.45127   47.12148    -1.58   0.114    -166.8521    17.94953
     g6_tavg |  -33.10528   44.19828    -0.75   0.454     -119.774    53.56339
     g7_tavg |   .4320024   42.99219     0.01   0.992    -83.87164    84.73565
     g8_tavg |   79.01014    47.1729     1.67   0.094     -13.4915    171.5118
     g9_tavg |   201.4064   50.52632     3.99   0.000      102.329    300.4837
     g1_wdir |  -274.5404   126.0237    -2.18   0.029    -521.6609    -27.4198
     g2_wdir |  -259.7535   127.1237    -2.04   0.041    -509.0312    -10.4759
     g3_wdir |  -418.6197   130.9358    -3.20   0.001    -675.3725    -161.867
     g4_wdir |   -445.263   127.8703    -3.48   0.001    -696.0046   -194.5214
     g5_wdir |  -389.4558   127.6333    -3.05   0.002    -639.7328   -139.1788
     g6_wdir |  -469.7882   138.2669    -3.40   0.001    -740.9166   -198.6598
    hms_deep |  -76.25929   18.86572    -4.04   0.000    -113.2532   -39.26538
       _cons |   32746.98   11167.63     2.93   0.003     10848.31    54645.66
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
          g_ppt |        10           0          10     |
         g_wspd |        10           1           9     |
 fe_countyqtroy |     11348       11348           0    *|
        fe_styr |       631           1         630    ?|
--------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", table(employment, `tbl
> _settings_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp g*_tavg g*_wdir hms_deep [aw=seer_pop] , a(g_ppt g_wspd fe_countyqtroy fe_styr) vce(clust
> er countyfip fe_stqtros)
(MWFE estimator converged in 15 iterations)

HDFE Linear regression                            Number of obs   =    147,352
Absorbing 4 HDFE groups                           F(  16,   2547) =      12.14
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2536
                                                  Adj R-squared   =     0.1873
Number of clusters (countyfip) =      2,837       Within R-sq.    =     0.0066
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 10843.3507

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     g1_tavg |  -61.50484   42.89498    -1.43   0.152    -145.6174    22.60774
     g2_tavg |  -91.43211   43.18815    -2.12   0.034    -176.1196   -6.744655
     g3_tavg |  -104.2987   30.19787    -3.45   0.001    -163.5136   -45.08382
     g4_tavg |  -46.59287   26.98096    -1.73   0.084    -99.49973    6.313987
     g5_tavg |  -48.41562   22.15833    -2.18   0.029    -91.86579   -4.965447
     g6_tavg |  -51.86442   18.61042    -2.79   0.005    -88.35752   -15.37133
     g7_tavg |    -38.515   20.11361    -1.91   0.056    -77.95569    .9256978
     g8_tavg |  -5.502354   20.45166    -0.27   0.788    -45.60592    34.60121
     g9_tavg |  -15.71719   24.30676    -0.65   0.518    -63.38021    31.94583
     g1_wdir |  -881.4541   83.07277   -10.61   0.000    -1044.351    -718.557
     g2_wdir |  -894.2912   83.28614   -10.74   0.000    -1057.607   -730.9757
     g3_wdir |  -909.0897    80.0763   -11.35   0.000    -1066.111   -752.0685
     g4_wdir |  -929.5637    82.6485   -11.25   0.000    -1091.629   -767.4986
     g5_wdir |  -916.7327   89.60645   -10.23   0.000    -1092.442   -741.0237
     g6_wdir |  -922.8525   84.39376   -10.94   0.000     -1088.34   -757.3651
    hms_deep |  -32.82859   9.589606    -3.42   0.001    -51.63281   -14.02437
       _cons |   83830.06   7375.435    11.37   0.000      69367.6    98292.52
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
          g_ppt |        10           0          10     |
         g_wspd |        10           1           9     |
 fe_countyqtroy |     11348       11348           0    *|
        fe_styr |       637           1         636    ?|
--------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", table(lfp, `tbl_settin
> gs_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(0 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.043***   -4.950***   -76.3***   -32.8*** |
  2. | hms_deep_stderr    (0.006)     (0.791)     (18.9)      (9.6) |
  3. |                                                              |
  4. |    observations     75,181     146,268    146,268    147,352 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b1.dta saved
.                 restore 
.                 
.         
.         }

.         
.         ** panel b, second row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(c.g1_wdir##i.statefip c.g2_wdir##i.statefip c.g3_wdir##i.statefip c.g4_wd
> ir##i.statefip c.g5_wdir##i.statefip c.g6_wdir##i.statefip fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 177 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 8 HDFE groups                           F(   1,   1685) =      55.69
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.8085
                                                  Adj R-squared   =     0.7869
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0287
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.2794

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0484104   .0064872     7.46   0.000     .0356866    .0611343
       _cons |   9.225685   .0370896   248.74   0.000     9.152938    9.298431
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------------+
        Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------------+---------------------------------------|
           statefip |        49           0          49     |
 statefip#c.g1_wdir |        49           0          49    ?|
           statefip |        49          49           0     |
 statefip#c.g2_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g3_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g4_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g5_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g6_wdir |        49           0          49    ?|
     fe_countyqtroy |      6677        6677           0    *|
            fe_styr |       637          49         588    ?|
------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", table(pm25, `tbl_setti
> ngs_log') replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(c.g1_wdir##i.statefip c.g2_wdir##i.statefip c.g3_wdir##i.stat
> efip c.g4_wdir##i.statefip c.g5_wdir##i.statefip c.g6_wdir##i.statefip fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 151 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 8 HDFE groups                           F(   1,   2519) =      40.03
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2539
                                                  Adj R-squared   =     0.1861
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0033
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   407.1251

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.018547   .7932366    -6.33   0.000     -6.57401   -3.463085
       _cons |   165.3723    4.53921    36.43   0.000     156.4714    174.2733
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------------+
        Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------------+---------------------------------------|
           statefip |        49           0          49     |
 statefip#c.g1_wdir |        49           0          49    ?|
           statefip |        49          49           0     |
 statefip#c.g2_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g3_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g4_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g5_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g6_wdir |        49           0          49    ?|
     fe_countyqtroy |     12424       12424           0    *|
            fe_styr |       631          49         582    ?|
------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", table(earnings, `tbl_s
> ettings_log') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(c.g1_wdir##i.statefip c.g2_wdir##i.statefip c.g3_wdi
> r##i.statefip c.g4_wdir##i.statefip c.g5_wdir##i.statefip c.g6_wdir##i.statefip fe_countyqtroy fe_styr) vce(cluster countyfip fe
> _stqtros)
(MWFE estimator converged in 151 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 8 HDFE groups                           F(   1,   2519) =       9.79
Statistics robust to heteroskedasticity           Prob > F        =     0.0018
                                                  R-squared       =     0.5378
                                                  Adj R-squared   =     0.4958
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0002
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 17874.4944

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -59.61089   19.04995    -3.13   0.002    -96.96606   -22.25572
       _cons |  -1496.304   117.1129   -12.78   0.000    -1725.952   -1266.657
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------------+
        Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------------+---------------------------------------|
           statefip |        49           0          49     |
 statefip#c.g1_wdir |        49           0          49    ?|
           statefip |        49          49           0     |
 statefip#c.g2_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g3_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g4_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g5_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g6_wdir |        49           0          49    ?|
     fe_countyqtroy |     12424       12424           0    *|
            fe_styr |       631          49         582    ?|
------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", table(employment, `tbl
> _settings_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(c.g1_wdir##i.statefip c.g2_wdir##i.statefip c.g3_wdir##i.state
> fip c.g4_wdir##i.statefip c.g5_wdir##i.statefip c.g6_wdir##i.statefip fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 158 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 8 HDFE groups                           F(   1,   2547) =      11.80
Statistics robust to heteroskedasticity           Prob > F        =     0.0006
                                                  R-squared       =     0.2523
                                                  Adj R-squared   =     0.1849
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0002
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11074.4326

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   -32.0207   9.321097    -3.44   0.001     -50.2984     -13.743
       _cons |  -1975.429   60.41714   -32.70   0.000    -2093.901   -1856.958
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------------+
        Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------------+---------------------------------------|
           statefip |        49           0          49     |
 statefip#c.g1_wdir |        49           0          49    ?|
           statefip |        49          49           0     |
 statefip#c.g2_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g3_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g4_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g5_wdir |        49           0          49    ?|
           statefip |        49          49           0    ?|
 statefip#c.g6_wdir |        49           0          49    ?|
     fe_countyqtroy |     12424       12424           0    *|
            fe_styr |       637          49         588    ?|
------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", table(lfp, `tbl_settin
> gs_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(0 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.048***   -5.019***   -59.6***   -32.0*** |
  2. | hms_deep_stderr    (0.006)     (0.793)     (19.0)      (9.3) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b2.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel b, third row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(c.g1_wdir##i.countyfip c.g2_wdir##i.countyfip c.g3_wdir##i.countyfip c.g4
> _wdir##i.countyfip c.g5_wdir##i.countyfip c.g6_wdir##i.countyfip fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 554 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 8 HDFE groups                           F(   1,   1685) =      49.90
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.8325
                                                  Adj R-squared   =     0.7358
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0300
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.4246

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0477392   .0067581     7.06   0.000     .0344841    .0609943
       _cons |   9.228943   .0373996   246.77   0.000     9.155588    9.302297
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------+
         Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------+---------------------------------------|
           countyfip |      1686        1686           0    *|
 countyfip#c.g1_wdir |      1686           2        1684    ?|
           countyfip |      1686        1686           0    *|
 countyfip#c.g2_wdir |      1686           1        1685    ?|
           countyfip |      1686        1686           0    *|
 countyfip#c.g3_wdir |      1686           0        1686    ?|
           countyfip |      1686        1686           0    *|
 countyfip#c.g4_wdir |      1686           4        1682     |
           countyfip |      1686        1686           0    *|
 countyfip#c.g5_wdir |      1686           2        1684    ?|
           countyfip |      1686        1686           0    *|
 countyfip#c.g6_wdir |      1686           1        1685    ?|
      fe_countyqtroy |      6677        6677           0    *|
             fe_styr |       637           0         637    ?|
-------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", table(pm25, `tbl_setti
> ngs_log') replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(c.g1_wdir##i.countyfip c.g2_wdir##i.countyfip c.g3_wdir##i.co
> untyfip c.g4_wdir##i.countyfip c.g5_wdir##i.countyfip c.g6_wdir##i.countyfip fe_countyqtroy fe_styr) vce(cluster countyfip fe_st
> qtros)
(MWFE estimator converged in 244 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 8 HDFE groups                           F(   1,   2519) =      37.79
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.3195
                                                  Adj R-squared   =     0.0086
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0038
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   449.3441

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   -5.31182   .8641274    -6.15   0.000    -7.006292   -3.617347
       _cons |   166.8554   4.787212    34.85   0.000     157.4682    176.2427
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------+
         Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------+---------------------------------------|
           countyfip |      3106        3106           0    *|
 countyfip#c.g1_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g2_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g3_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g4_wdir |      3106           6        3100     |
           countyfip |      3106        3106           0    *|
 countyfip#c.g5_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g6_wdir |      3106           6        3100    ?|
      fe_countyqtroy |     12424       12424           0    *|
             fe_styr |       631           0         631    ?|
-------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", table(earnings, `tbl_s
> ettings_log') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(c.g1_wdir##i.countyfip c.g2_wdir##i.countyfip c.g3_w
> dir##i.countyfip c.g4_wdir##i.countyfip c.g5_wdir##i.countyfip c.g6_wdir##i.countyfip fe_countyqtroy fe_styr) vce(cluster county
> fip fe_stqtros)
(MWFE estimator converged in 246 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 8 HDFE groups                           F(   1,   2519) =       9.68
Statistics robust to heteroskedasticity           Prob > F        =     0.0019
                                                  R-squared       =     0.5938
                                                  Adj R-squared   =     0.4082
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 19364.6574

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -63.10821   20.28072    -3.11   0.002    -102.8768   -23.33961
       _cons |   -1478.63   125.9113   -11.74   0.000    -1725.531    -1231.73
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------+
         Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------+---------------------------------------|
           countyfip |      3106        3106           0    *|
 countyfip#c.g1_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g2_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g3_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g4_wdir |      3106           6        3100     |
           countyfip |      3106        3106           0    *|
 countyfip#c.g5_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g6_wdir |      3106           6        3100    ?|
      fe_countyqtroy |     12424       12424           0    *|
             fe_styr |       631           0         631    ?|
-------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", table(employment, `tbl
> _settings_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(c.g1_wdir##i.countyfip c.g2_wdir##i.countyfip c.g3_wdir##i.cou
> ntyfip c.g4_wdir##i.countyfip c.g5_wdir##i.countyfip c.g6_wdir##i.countyfip fe_countyqtroy fe_styr) vce(cluster countyfip fe_stq
> tros)
(MWFE estimator converged in 246 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 8 HDFE groups                           F(   1,   2547) =      10.48
Statistics robust to heteroskedasticity           Prob > F        =     0.0012
                                                  R-squared       =     0.3662
                                                  Adj R-squared   =     0.0795
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0002
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11768.4513

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -31.99774    9.88289    -3.24   0.001    -51.37705   -12.61842
       _cons |  -1975.545   62.16173   -31.78   0.000    -2097.438   -1853.653
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------+
         Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------+---------------------------------------|
           countyfip |      3106        3106           0    *|
 countyfip#c.g1_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g2_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g3_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g4_wdir |      3106           6        3100     |
           countyfip |      3106        3106           0    *|
 countyfip#c.g5_wdir |      3106           6        3100    ?|
           countyfip |      3106        3106           0    *|
 countyfip#c.g6_wdir |      3106           6        3100    ?|
      fe_countyqtroy |     12424       12424           0    *|
             fe_styr |       637           0         637    ?|
-------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave hms_deep using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", table(lfp, `tbl_settin
> gs_pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(0 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.048***   -5.312***   -63.1***   -32.0*** |
  2. | hms_deep_stderr    (0.007)     (0.864)     (20.3)      (9.9) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_b3.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel c, first row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      71.70
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0369
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.3481

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0065709     8.47   0.000      .042752    .0685282
       _cons |   9.190593    .040427   227.34   0.000     9.111301    9.269886
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      45.21
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0038
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   409.6540

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   .7758234    -6.72   0.000    -6.737943   -3.695309
       _cons |    166.374   4.714186    35.29   0.000       157.13    175.6181
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stq
> tros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      13.20
Statistics robust to heteroskedasticity           Prob > F        =     0.0003
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0005
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 18020.2635

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   21.90873    -3.63   0.000    -122.5488   -36.62692
       _cons |  -1395.351   134.9056   -10.34   0.000    -1659.888   -1130.813
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   2547) =      17.57
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11151.6317

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033   9.231539    -4.19   0.000    -56.80242   -20.59825
       _cons |  -1941.697   64.50074   -30.10   0.000    -2068.176   -1815.217
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217***   -79.6***   -38.7*** |
  2. | hms_deep_stderr    (0.007)     (0.776)     (21.9)      (9.2) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c1.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel c, second row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_divisionyr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      62.22
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7611
                                                  Adj R-squared   =     0.7374
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0319
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.4203

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0532975   .0067567     7.89   0.000     .0400451      .06655
       _cons |   9.201964   .0426411   215.80   0.000     9.118329    9.285599
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
  fe_divisionyr |       117           0         117     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_divisionyr) vce(cluster countyfip fe_stqtro
> s)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      50.62
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2271
                                                  Adj R-squared   =     0.1615
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0042
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   413.2324

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.399253    .758842    -7.12   0.000    -6.887271   -3.911235
       _cons |   167.2976   4.738522    35.31   0.000     158.0058    176.5894
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
  fe_divisionyr |       117           0         117     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_divisionyr) vce(cluster countyfip 
> fe_stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      14.78
Statistics robust to heteroskedasticity           Prob > F        =     0.0001
                                                  R-squared       =     0.5013
                                                  Adj R-squared   =     0.4589
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0007
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 18516.3232

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -100.7105    26.1931    -3.84   0.000    -152.0727   -49.34824
       _cons |  -1288.608    173.368    -7.43   0.000    -1628.566   -948.6496
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
  fe_divisionyr |       117           0         117     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_divisionyr) vce(cluster countyfip fe_stqtros
> )
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   2547) =       2.53
Statistics robust to heteroskedasticity           Prob > F        =     0.1117
                                                  R-squared       =     0.1579
                                                  Adj R-squared   =     0.0870
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11720.2536

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -27.86031   17.50958    -1.59   0.112    -62.19477     6.47414
       _cons |   -1996.44   122.4209   -16.31   0.000    -2236.494   -1756.385
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
  fe_divisionyr |       117           0         117     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employm~t       lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.053***   -5.399***   -100.7***     -27.9 |
  2. | hms_deep_stderr    (0.007)     (0.759)      (26.2)    (17.5) |
  3. |                                                              |
  4. |    observations     75,207     160,346     160,346   161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c2.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel c, third row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_regionyr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      61.61
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7527
                                                  Adj R-squared   =     0.7284
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0332
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.4444

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0542433   .0069108     7.85   0.000     .0406886     .067798
       _cons |   9.197373    .042682   215.49   0.000     9.113658    9.281089
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
    fe_regionyr |        52           0          52     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_regionyr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      45.50
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2215
                                                  Adj R-squared   =     0.1559
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0041
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   414.6244

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.212464   .7727812    -6.75   0.000    -6.727815   -3.697112
       _cons |    166.353   4.808422    34.60   0.000     156.9241    175.7818
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
    fe_regionyr |        52           0          52     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_regionyr) vce(cluster countyfip fe
> _stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      13.59
Statistics robust to heteroskedasticity           Prob > F        =     0.0002
                                                  R-squared       =     0.4888
                                                  Adj R-squared   =     0.4456
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0009
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 18742.5594

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -109.2334   29.63508    -3.69   0.000     -167.345   -51.12181
       _cons |  -1245.537   201.7776    -6.17   0.000    -1641.204   -849.8704
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
    fe_regionyr |        52           0          52     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_regionyr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   2547) =       1.61
Statistics robust to heteroskedasticity           Prob > F        =     0.2052
                                                  R-squared       =     0.1315
                                                  Adj R-squared   =     0.0588
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11899.7824

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -26.02048    20.5328    -1.27   0.205    -66.28317    14.24221
       _cons |  -2005.731   143.8208   -13.95   0.000    -2287.748   -1723.713
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
    fe_regionyr |        52           0          52     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employm~t       lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.054***   -5.212***   -109.2***     -26.0 |
  2. | hms_deep_stderr    (0.007)     (0.773)      (29.6)    (20.5) |
  3. |                                                              |
  4. |    observations     75,207     160,346     160,346   161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c3.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel c, fourth row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy rfrnc_yr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      47.66
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7309
                                                  Adj R-squared   =     0.7046
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0370
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.5064

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0568862     .00824     6.90   0.000     .0407244    .0730479
       _cons |   9.184545   .0470525   195.20   0.000     9.092258    9.276833
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
       rfrnc_yr |        13           0          13     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy rfrnc_yr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      30.79
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2126
                                                  Adj R-squared   =     0.1464
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0032
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   416.9361

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -4.421408   .7968562    -5.55   0.000    -5.983969   -2.858848
       _cons |   162.3525    5.18383    31.32   0.000     152.1875    172.5175
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
       rfrnc_yr |        13           0          13     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy rfrnc_yr) vce(cluster countyfip fe_st
> qtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   2519) =      15.99
Statistics robust to heteroskedasticity           Prob > F        =     0.0001
                                                  R-squared       =     0.4736
                                                  Adj R-squared   =     0.4293
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0014
Number of clusters (fe_stqtros) =      2,520      Root MSE        = 19017.0990

               (Std. Err. adjusted for 2,520 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -132.9057   33.23192    -4.00   0.000    -198.0704   -67.74102
       _cons |   -1125.91   228.7624    -4.92   0.000    -1574.491   -677.3282
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
       rfrnc_yr |        13           0          13     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy rfrnc_yr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   2547) =       0.37
Statistics robust to heteroskedasticity           Prob > F        =     0.5437
                                                  R-squared       =     0.1096
                                                  Adj R-squared   =     0.0353
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0000
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 12047.7629

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   13.54641   22.30422     0.61   0.544    -30.18983    57.28266
       _cons |  -2205.547   159.6033   -13.82   0.000    -2518.513   -1892.582
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
       rfrnc_yr |        13           0          13     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employm~t       lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.057***   -4.421***   -132.9***      13.5 |
  2. | hms_deep_stderr    (0.008)     (0.797)      (33.2)    (22.3) |
  3. |                                                              |
  4. |    observations     75,207     160,346     160,346   161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_c4.dta saved
.                 restore 
.                 
.         }

.         
.         
.         use "$Rep_smokelabor/1_build/regdata/county_year.dta", clear 

.         
.         local tbl_settings_log format(%6.3f) parentheses(stderr) asterisk()

.         local tbl_settings_pmil format(%6.1f) parentheses(stderr) asterisk()

.         
.         ** panel d, first row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(countyfip fe_styr) vce(cluster countyfip fe_styr)
(dropped 54 singleton observations)
(MWFE estimator converged in 8 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     19,011
Absorbing 2 HDFE groups                           F(   1,    618) =       7.15
Statistics robust to heteroskedasticity           Prob > F        =     0.0077
                                                  R-squared       =     0.9218
                                                  Adj R-squared   =     0.9110
Number of clusters (countyfip) =      1,688       Within R-sq.    =     0.0372
Number of clusters (fe_styr) =        619         Root MSE        =     0.6676

                    (Std. Err. adjusted for 619 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0359903   .0134595     2.67   0.008     .0095585    .0624222
       _cons |   8.758214   .2615288    33.49   0.000     8.244621    9.271807
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      1688        1688           0    *|
     fe_styr |       619         619           0    *|
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta saved
.                 sleep 3000 
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(countyfip fe_styr) vce(cluster countyfip fe_styr)
(dropped 13 singleton observations)
(MWFE estimator converged in 4 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     40,152
Absorbing 2 HDFE groups                           F(   1,    617) =       0.04
Statistics robust to heteroskedasticity           Prob > F        =     0.8406
                                                  R-squared       =     0.4701
                                                  Adj R-squared   =     0.4159
Number of clusters (countyfip) =      3,105       Within R-sq.    =     0.0000
Number of clusters (fe_styr) =        618         Root MSE        =  1038.9646

                    (Std. Err. adjusted for 618 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   1.083213    5.38317     0.20   0.841    -9.488345    11.65477
       _cons |   524.3953   106.9105     4.90   0.000     314.4427    734.3478
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3105        3105           0    *|
     fe_styr |       618         618           0    *|
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta saved
.                 sleep 3000 
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(countyfip fe_styr) vce(cluster countyfip fe_styr)
(dropped 13 singleton observations)
(MWFE estimator converged in 4 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     40,152
Absorbing 2 HDFE groups                           F(   1,    617) =       0.01
Statistics robust to heteroskedasticity           Prob > F        =     0.9353
                                                  R-squared       =     0.6515
                                                  Adj R-squared   =     0.6159
Number of clusters (countyfip) =      3,105       Within R-sq.    =     0.0000
Number of clusters (fe_styr) =        618         Root MSE        = 14034.2658

                    (Std. Err. adjusted for 618 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -3.538481   43.56818    -0.08   0.935    -89.09838    82.02142
       _cons |  -1723.118   877.5121    -1.96   0.050    -3446.391    .1541036
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3105        3105           0    *|
     fe_styr |       618         618           0    *|
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta saved
.                 sleep 3000 
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(countyfip fe_styr) vce(cluster countyfip fe_styr)
(dropped 13 singleton observations)
(MWFE estimator converged in 5 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     40,358
Absorbing 2 HDFE groups                           F(   1,    623) =       0.75
Statistics robust to heteroskedasticity           Prob > F        =     0.3871
                                                  R-squared       =     0.2824
                                                  Adj R-squared   =     0.2094
Number of clusters (countyfip) =      3,105       Within R-sq.    =     0.0001
Number of clusters (fe_styr) =        624         Root MSE        = 10028.2427

                    (Std. Err. adjusted for 624 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -30.22419   34.91959    -0.87   0.387    -98.79856    38.35017
       _cons |  -1535.076    708.089    -2.17   0.031    -2925.606   -144.5455
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3105        3105           0    *|
     fe_styr |       624         624           0    *|
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta saved
.                 sleep 3000 
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +-----------------------------------------------------------+
     |             var       pm25   earnings   employ~t      lfp |
     |-----------------------------------------------------------|
  1. |   hms_deep_coef   0.036***      1.083       -3.5    -30.2 |
  2. | hms_deep_stderr    (0.013)    (5.383)     (43.6)   (34.9) |
  3. |                                                           |
  4. |    observations     19,011     40,152     40,152   40,358 |
     +-----------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d1.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel d, second row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(countyfip fe_divisionyr) vce(cluster countyfip fe_styr)
(dropped 36 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =     19,029
Absorbing 2 HDFE groups                           F(   1,    636) =       3.98
Statistics robust to heteroskedasticity           Prob > F        =     0.0465
                                                  R-squared       =     0.8811
                                                  Adj R-squared   =     0.8686
Number of clusters (countyfip) =      1,689       Within R-sq.    =     0.0111
Number of clusters (fe_styr) =        637         Root MSE        =     0.8105

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0157982   .0079191     1.99   0.046     .0002475     .031349
       _cons |   9.151037   .1550589    59.02   0.000     8.846547    9.455526
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------+
   Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------+---------------------------------------|
     countyfip |      1689        1689           0    *|
 fe_divisionyr |       117           0         117     |
-------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta saved
.                 sleep 3000 
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(countyfip fe_divisionyr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 5 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       4.25
Statistics robust to heteroskedasticity           Prob > F        =     0.0397
                                                  R-squared       =     0.3982
                                                  Adj R-squared   =     0.3457
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0009
Number of clusters (fe_styr) =        631         Root MSE        =  1099.9189

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.902132   2.862881    -2.06   0.040    -11.52408   -.2801878
       _cons |   667.5805   55.65487    12.00   0.000      558.289     776.872
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------+
   Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------+---------------------------------------|
     countyfip |      3106        3106           0    *|
 fe_divisionyr |       117           0         117     |
-------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta saved
.                 sleep 3000 
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(countyfip fe_divisionyr) vce(cluster countyfip fe_st
> yr)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       5.78
Statistics robust to heteroskedasticity           Prob > F        =     0.0165
                                                  R-squared       =     0.6176
                                                  Adj R-squared   =     0.5842
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0009
Number of clusters (fe_styr) =        631         Root MSE        = 14597.0195

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -80.00589   33.27767    -2.40   0.016    -145.3545   -14.65731
       _cons |  -179.7103   684.0972    -0.26   0.793    -1523.097    1163.676
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------+
   Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------+---------------------------------------|
     countyfip |      3106        3106           0    *|
 fe_divisionyr |       117           0         117     |
-------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta saved
.                 sleep 3000 
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(countyfip fe_divisionyr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =     40,371
Absorbing 2 HDFE groups                           F(   1,    636) =       0.37
Statistics robust to heteroskedasticity           Prob > F        =     0.5453
                                                  R-squared       =     0.1854
                                                  Adj R-squared   =     0.1147
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_styr) =        637         Root MSE        = 10604.9145

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   15.77109   26.05907     0.61   0.545    -35.40113    66.94331
       _cons |  -2454.415   539.8351    -4.55   0.000     -3514.49   -1394.341
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------+
   Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------+---------------------------------------|
     countyfip |      3106        3106           0    *|
 fe_divisionyr |       117           0         117     |
-------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta saved
.                 sleep 3000 
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +----------------------------------------------------------+
     |             var      pm25   earnings   employ~t      lfp |
     |----------------------------------------------------------|
  1. |   hms_deep_coef   0.016**   -5.902**    -80.0**     15.8 |
  2. | hms_deep_stderr   (0.008)    (2.863)     (33.3)   (26.1) |
  3. |                                                          |
  4. |    observations    19,029     40,165     40,165   40,371 |
     +----------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d2.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel d, third row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(countyfip fe_regionyr) vce(cluster countyfip fe_styr)
(dropped 36 singleton observations)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     19,029
Absorbing 2 HDFE groups                           F(   1,    636) =       7.18
Statistics robust to heteroskedasticity           Prob > F        =     0.0076
                                                  R-squared       =     0.8683
                                                  Adj R-squared   =     0.8550
Number of clusters (countyfip) =      1,689       Within R-sq.    =     0.0179
Number of clusters (fe_styr) =        637         Root MSE        =     0.8515

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0169857   .0063403     2.68   0.008     .0045352    .0294362
       _cons |   9.127985   .1214276    75.17   0.000     8.889538    9.366433
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      1689        1689           0    *|
 fe_regionyr |        52           0          52     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta saved
.                 sleep 3000 
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(countyfip fe_regionyr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       3.90
Statistics robust to heteroskedasticity           Prob > F        =     0.0487
                                                  R-squared       =     0.3843
                                                  Adj R-squared   =     0.3318
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0007
Number of clusters (fe_styr) =        631         Root MSE        =  1111.5262

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -4.139335   2.096126    -1.97   0.049    -8.255574   -.0230957
       _cons |   631.8945   40.05135    15.78   0.000     553.2443    710.5448
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
 fe_regionyr |        52           0          52     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta saved
.                 sleep 3000 
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(countyfip fe_regionyr) vce(cluster countyfip fe_styr
> )
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       4.21
Statistics robust to heteroskedasticity           Prob > F        =     0.0406
                                                  R-squared       =     0.6021
                                                  Adj R-squared   =     0.5682
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0014
Number of clusters (fe_styr) =        631         Root MSE        = 14875.0159

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.09945   38.54319    -2.05   0.041    -154.7881    -3.41077
       _cons |  -198.0469   784.0421    -0.25   0.801    -1737.699    1341.605
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
 fe_regionyr |        52           0          52     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta saved
.                 sleep 3000 
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(countyfip fe_regionyr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =     40,371
Absorbing 2 HDFE groups                           F(   1,    636) =       0.20
Statistics robust to heteroskedasticity           Prob > F        =     0.6547
                                                  R-squared       =     0.1542
                                                  Adj R-squared   =     0.0825
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_styr) =        637         Root MSE        = 10796.4830

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   11.16468   24.95375     0.45   0.655    -37.83702    60.16639
       _cons |  -2361.361   517.1545    -4.57   0.000    -3376.898   -1345.825
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
 fe_regionyr |        52           0          52     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta saved
.                 sleep 3000 
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +-----------------------------------------------------------+
     |             var       pm25   earnings   employ~t      lfp |
     |-----------------------------------------------------------|
  1. |   hms_deep_coef   0.017***   -4.139**    -79.1**     11.2 |
  2. | hms_deep_stderr    (0.006)    (2.096)     (38.5)   (25.0) |
  3. |                                                           |
  4. |    observations     19,029     40,165     40,165   40,371 |
     +-----------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d3.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel d, fourth row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(countyfip rfrnc_yr) vce(cluster countyfip fe_styr)
(dropped 36 singleton observations)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =     19,029
Absorbing 2 HDFE groups                           F(   1,    636) =      11.46
Statistics robust to heteroskedasticity           Prob > F        =     0.0008
                                                  R-squared       =     0.8350
                                                  Adj R-squared   =     0.8188
Number of clusters (countyfip) =      1,689       Within R-sq.    =     0.0325
Number of clusters (fe_styr) =        637         Root MSE        =     0.9518

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0187128   .0055269     3.39   0.001     .0078596    .0295659
       _cons |    9.09446   .1039902    87.45   0.000     8.890254    9.298666
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      1689        1689           0    *|
    rfrnc_yr |        13           0          13     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta saved
.                 sleep 3000 
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(countyfip rfrnc_yr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       0.00
Statistics robust to heteroskedasticity           Prob > F        =     0.9846
                                                  R-squared       =     0.3673
                                                  Adj R-squared   =     0.3141
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0000
Number of clusters (fe_styr) =        631         Root MSE        =  1126.1893

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |    .044859   2.317488     0.02   0.985    -4.506076    4.595794
       _cons |     547.19   47.06232    11.63   0.000      454.772     639.608
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
    rfrnc_yr |        13           0          13     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta saved
.                 sleep 3000 
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(countyfip rfrnc_yr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =     40,165
Absorbing 2 HDFE groups                           F(   1,    630) =       6.50
Statistics robust to heteroskedasticity           Prob > F        =     0.0110
                                                  R-squared       =     0.5836
                                                  Adj R-squared   =     0.5486
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0027
Number of clusters (fe_styr) =        631         Root MSE        = 15209.4636

                    (Std. Err. adjusted for 631 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -82.82505    32.4744    -2.55   0.011    -146.5962   -19.05389
       _cons |  -122.6807   665.8656    -0.18   0.854    -1430.265    1184.904
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
    rfrnc_yr |        13           0          13     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta saved
.                 sleep 3000 
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(countyfip rfrnc_yr) vce(cluster countyfip fe_styr)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =     40,371
Absorbing 2 HDFE groups                           F(   1,    636) =       6.57
Statistics robust to heteroskedasticity           Prob > F        =     0.0106
                                                  R-squared       =     0.1301
                                                  Adj R-squared   =     0.0572
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0020
Number of clusters (fe_styr) =        637         Root MSE        = 10943.8612

                    (Std. Err. adjusted for 637 clusters in countyfip fe_styr)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |     51.249   19.98809     2.56   0.011     11.99837    90.49963
       _cons |  -3171.105   455.4207    -6.96   0.000    -4065.415   -2276.795
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
   countyfip |      3106        3106           0    *|
    rfrnc_yr |        13           0          13     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta saved
.                 sleep 3000 
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +-----------------------------------------------------------+
     |             var       pm25   earnings   employ~t      lfp |
     |-----------------------------------------------------------|
  1. |   hms_deep_coef   0.019***      0.045    -82.8**   51.2** |
  2. | hms_deep_stderr    (0.006)    (2.317)     (32.5)   (20.0) |
  3. |                                                           |
  4. |    observations     19,029     40,165     40,165   40,371 |
     +-----------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_d4.dta saved
.                 restore 
.                 
.         }

.         
.         use "$Rep_smokelabor/1_build/regdata/county_quarter.dta", clear 

.         
.         local tbl_settings_log format(%6.3f) parentheses(stderr) asterisk()

.         local tbl_settings_pmil format(%6.1f) parentheses(stderr) asterisk()

.         
.         ** panel e, first row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr i.fe_countyqtroy##c.rfrnc_yr) vce(cluster countyfi
> p fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 9 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 3 HDFE groups                           F(   1,   1685) =     116.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.8307
                                                  Adj R-squared   =     0.7665
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0463
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.3392

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0595432    .005511    10.80   0.000     .0487341    .0703522
       _cons |   9.171649   .0374016   245.22   0.000      9.09829    9.245007
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------------+
               Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------------+---------------------------------------|
            fe_countyqtroy |      6677        6677           0    *|
                   fe_styr |       637           0         637     |
            fe_countyqtroy |      6677        6677           0    *|
 fe_countyqtroy#c.rfrnc_yr |      6677           0        6677    ?|
-------------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta saved
.                 
.                 * income 
.                 reghdfe pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr i.fe_countyqtroy##c.rfrnc_yr) vce(cluste
> r countyfip fe_stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,402
Absorbing 3 HDFE groups                           F(   1,   2523) =      16.93
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9913
                                                  Adj R-squared   =     0.9887
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0010
Number of clusters (fe_stqtros) =      2,524      Root MSE        =   391.2326

               (Std. Err. adjusted for 2,524 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pc_qwi_pay~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -2.453764   .5962778    -4.12   0.000    -3.623008    -1.28452
       _cons |   5374.214   3.290294  1633.35   0.000     5367.762    5380.666
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------------+
               Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------------+---------------------------------------|
            fe_countyqtroy |     12424       12424           0    *|
                   fe_styr |       632           0         632     |
            fe_countyqtroy |     12424       12424           0    *|
 fe_countyqtroy#c.rfrnc_yr |     12424           0       12424    ?|
-------------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta saved
.                 
.                 * emp 
.                 reghdfe pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr i.fe_countyqtroy##c.rfrnc_yr) v
> ce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,402
Absorbing 3 HDFE groups                           F(   1,   2523) =      14.87
Statistics robust to heteroskedasticity           Prob > F        =     0.0001
                                                  R-squared       =     0.9937
                                                  Adj R-squared   =     0.9918
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,524      Root MSE        = 19059.7644

               (Std. Err. adjusted for 2,524 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pmil_qwi_e~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   -42.6865    11.0693    -3.86   0.000    -64.39233   -20.98066
       _cons |   626047.2    76.4574  8188.18   0.000     625897.3    626197.1
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------------+
               Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------------+---------------------------------------|
            fe_countyqtroy |     12424       12424           0    *|
                   fe_styr |       632           0         632     |
            fe_countyqtroy |     12424       12424           0    *|
 fe_countyqtroy#c.rfrnc_yr |     12424           0       12424    ?|
-------------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta saved
.                 
.                 * lfp
.                 reghdfe pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr i.fe_countyqtroy##c.rfrnc_yr) vce(cluster
>  countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =    161,512
Absorbing 3 HDFE groups                           F(   1,   2547) =      11.15
Statistics robust to heteroskedasticity           Prob > F        =     0.0009
                                                  R-squared       =     0.9710
                                                  Adj R-squared   =     0.9621
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 11867.8611

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pmil_lau_lfp |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -24.73577   7.406931    -3.34   0.001    -39.25999   -10.21155
       _cons |   625555.7   26.84812  2.3e+04   0.000       625503    625608.3
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-------------------------------------------------------------------+
               Absorbed FE | Categories  - Redundant  = Num. Coefs |
---------------------------+---------------------------------------|
            fe_countyqtroy |     12424       12424           0    *|
                   fe_styr |       637           0         637     |
            fe_countyqtroy |     12424       12424           0    *|
 fe_countyqtroy#c.rfrnc_yr |     12424           0       12424    ?|
-------------------------------------------------------------------+
? = number of redundant parameters may be higher
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.060***   -2.454***   -42.7***   -24.7*** |
  2. | hms_deep_stderr    (0.006)     (0.596)     (11.1)      (7.4) |
  3. |                                                              |
  4. |    observations     75,207     160,402    160,402    161,512 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e1.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel e, second row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =      71.70
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0369
Number of clusters (fe_stqtros) =      2,548      Root MSE        =     1.3481

               (Std. Err. adjusted for 1,686 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0065709     8.47   0.000      .042752    .0685282
       _cons |   9.190593    .040427   227.34   0.000     9.111301    9.269886
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta saved
.                 
.                 * income 
.                 reghdfe pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 5 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =    160,402
Absorbing 2 HDFE groups                           F(   1,   2523) =       2.68
Statistics robust to heteroskedasticity           Prob > F        =     0.1015
                                                  R-squared       =     0.9728
                                                  Adj R-squared   =     0.9704
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,524      Root MSE        =   631.9194

               (Std. Err. adjusted for 2,524 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pc_qwi_pay~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -1.581239   .9652214    -1.64   0.102    -3.473946    .3114679
       _cons |   5369.802   4.452128  1206.12   0.000     5361.071    5378.532
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       632           0         632     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta saved
.                 
.                 * emp 
.                 reghdfe pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtr
> os)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,402
Absorbing 2 HDFE groups                           F(   1,   2523) =       1.57
Statistics robust to heteroskedasticity           Prob > F        =     0.2105
                                                  R-squared       =     0.9869
                                                  Adj R-squared   =     0.9858
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0000
Number of clusters (fe_stqtros) =      2,524      Root MSE        = 25088.7676

               (Std. Err. adjusted for 2,524 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pmil_qwi_e~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   10.76746   8.597295     1.25   0.211    -6.091021    27.62593
       _cons |   625777.1   16.33441  3.8e+04   0.000       625745    625809.1
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       632           0         632     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta saved
.                 
.                 * lfp
.                 reghdfe pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =    161,512
Absorbing 2 HDFE groups                           F(   1,   2547) =     151.53
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9344
                                                  Adj R-squared   =     0.9286
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0001
Number of clusters (fe_stqtros) =      2,548      Root MSE        = 16297.0683

               (Std. Err. adjusted for 2,548 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
pmil_lau_lfp |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -25.34189    2.05866   -12.31   0.000    -29.37871   -21.30507
       _cons |   625558.8   2.476568  2.5e+05   0.000     625553.9    625563.6
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +-------------------------------------------------------------+
     |             var       pm25   earnings   employ~t        lfp |
     |-------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***     -1.581       10.8   -25.3*** |
  2. | hms_deep_stderr    (0.007)    (0.965)      (8.6)      (2.1) |
  3. |                                                             |
  4. |    observations     75,207    160,402    160,402    161,512 |
     +-------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e2.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel e, third row 
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe d_pm25 d_hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(dropped 110 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     67,536
Absorbing 2 HDFE groups                           F(   1,   1652) =      99.16
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2478
                                                  Adj R-squared   =     0.1587
Number of clusters (countyfip) =      1,653       Within R-sq.    =     0.0467
Number of clusters (fe_stqtros) =      2,352      Root MSE        =     1.7617

               (Std. Err. adjusted for 1,653 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
      d_pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  d_hms_deep |   .0629251   .0063192     9.96   0.000     .0505307    .0753195
       _cons |  -.3448305    .037077    -9.30   0.000    -.4175534   -.2721077
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6566        6566           0    *|
        fe_styr |       588           0         588     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll d_hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    147,978
Absorbing 2 HDFE groups                           F(   1,   2327) =      31.39
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2373
                                                  Adj R-squared   =     0.1638
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0031
Number of clusters (fe_stqtros) =      2,328      Root MSE        =   414.3400

               (Std. Err. adjusted for 2,328 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  d_hms_deep |  -3.662924   .6537808    -5.60   0.000    -4.944978    -2.38087
       _cons |   135.7502   2.897548    46.85   0.000     130.0681    141.4322
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       583           0         583     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal d_hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_s
> tqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    147,978
Absorbing 2 HDFE groups                           F(   1,   2327) =      10.24
Statistics robust to heteroskedasticity           Prob > F        =     0.0014
                                                  R-squared       =     0.5447
                                                  Adj R-squared   =     0.5009
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (fe_stqtros) =      2,328      Root MSE        = 18040.8294

               (Std. Err. adjusted for 2,328 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  d_hms_deep |  -51.57996   16.11994    -3.20   0.001    -83.19091   -19.96901
       _cons |  -1686.242   129.3228   -13.04   0.000    -1939.842   -1432.642
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       583           0         583     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp d_hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_stqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    149,088
Absorbing 2 HDFE groups                           F(   1,   2351) =       8.85
Statistics robust to heteroskedasticity           Prob > F        =     0.0030
                                                  R-squared       =     0.2447
                                                  Adj R-squared   =     0.1725
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0002
Number of clusters (fe_stqtros) =      2,352      Root MSE        = 11244.9422

               (Std. Err. adjusted for 2,352 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  d_hms_deep |  -26.86678   9.030021    -2.98   0.003    -44.57441   -9.159145
       _cons |  -2197.002   60.80711   -36.13   0.000    -2316.243   -2077.761
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       588           0         588     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +----------------------------------------------------------------+
     |               var       pm25    earnings   employ~t        lfp |
     |----------------------------------------------------------------|
  1. |   d_hms_deep_coef   0.063***   -3.663***   -51.6***   -26.9*** |
  2. | d_hms_deep_stderr    (0.006)     (0.654)     (16.1)      (9.0) |
  3. |                                                                |
  4. |      observations     67,536     147,978    147,978    149,088 |
     +----------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_e3.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel f, first row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_divisionqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,    467) =      35.79
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0369
Number of clusters (fe_divisionqtros) =        468Root MSE        =     1.3481

           (Std. Err. adjusted for 468 clusters in countyfip fe_divisionqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0093011     5.98   0.000     .0373628    .0739174
       _cons |   9.190593   .0567468   161.96   0.000     9.079083    9.302104
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_divisionqtro
> s)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,    467) =      24.11
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0038
Number of clusters (fe_divisionqtros) =        468Root MSE        =   409.6540

           (Std. Err. adjusted for 468 clusters in countyfip fe_divisionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   1.062316    -4.91   0.000    -7.304138   -3.129114
       _cons |    166.374    6.85078    24.29   0.000     152.9119    179.8362
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_div
> isionqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,    467) =       4.46
Statistics robust to heteroskedasticity           Prob > F        =     0.0353
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0005
Number of clusters (fe_divisionqtros) =        468Root MSE        = 18020.2635

           (Std. Err. adjusted for 468 clusters in countyfip fe_divisionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   37.70266    -2.11   0.035    -153.6757   -5.500014
       _cons |  -1395.351   251.0412    -5.56   0.000    -1888.661   -902.0404
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_divisionqtros
> )
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,    467) =       8.76
Statistics robust to heteroskedasticity           Prob > F        =     0.0032
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (fe_divisionqtros) =        468Root MSE        = 11151.6317

           (Std. Err. adjusted for 468 clusters in countyfip fe_divisionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033   13.07371    -2.96   0.003    -64.39091   -13.00976
       _cons |  -1941.697   97.68327   -19.88   0.000     -2133.65   -1749.743
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217***    -79.6**   -38.7*** |
  2. | hms_deep_stderr    (0.009)     (1.062)     (37.7)     (13.1) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f1.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel f, second row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_regionqtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,    207) =      25.90
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0369
Number of clusters (fe_regionqtros) =        208  Root MSE        =     1.3481

             (Std. Err. adjusted for 208 clusters in countyfip fe_regionqtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0109319     5.09   0.000     .0340879    .0771923
       _cons |   9.190593   .0669942   137.18   0.000     9.058515    9.322672
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_regionqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,    207) =      13.26
Statistics robust to heteroskedasticity           Prob > F        =     0.0003
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0038
Number of clusters (fe_regionqtros) =        208  Root MSE        =   409.6540

             (Std. Err. adjusted for 208 clusters in countyfip fe_regionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   1.432459    -3.64   0.000    -8.040706   -2.392546
       _cons |    166.374   9.420699    17.66   0.000     147.8012    184.9468
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_reg
> ionqtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,    207) =       2.35
Statistics robust to heteroskedasticity           Prob > F        =     0.1270
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0005
Number of clusters (fe_regionqtros) =        208  Root MSE        = 18020.2635

             (Std. Err. adjusted for 208 clusters in countyfip fe_regionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   51.93703    -1.53   0.127    -181.9812    22.80548
       _cons |  -1395.351   350.0667    -3.99   0.000    -2085.504   -705.1976
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_regionqtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,    207) =       5.54
Statistics robust to heteroskedasticity           Prob > F        =     0.0196
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (fe_regionqtros) =        208  Root MSE        = 11151.6317

             (Std. Err. adjusted for 208 clusters in countyfip fe_regionqtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033   16.44808    -2.35   0.020    -71.12756   -6.273105
       _cons |  -1941.697   132.7314   -14.63   0.000    -2203.375   -1680.018
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +-------------------------------------------------------------+
     |             var       pm25    earnings   employ~t       lfp |
     |-------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217***      -79.6   -38.7** |
  2. | hms_deep_stderr    (0.011)     (1.432)     (51.9)    (16.4) |
  3. |                                                             |
  4. |    observations     75,207     160,346    160,346   161,498 |
     +-------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f2.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel f, third row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip rfrnc_qtros)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,     51) =      24.62
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
Number of clusters (countyfip) =      1,686       Within R-sq.    =     0.0369
Number of clusters (rfrnc_qtros) =         52     Root MSE        =     1.3481

                 (Std. Err. adjusted for 52 clusters in countyfip rfrnc_qtros)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0112142     4.96   0.000     .0331267    .0781535
       _cons |   9.190593   .0860135   106.85   0.000     9.017914    9.363273
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip rfrnc_qtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,     51) =       6.31
Statistics robust to heteroskedasticity           Prob > F        =     0.0152
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0038
Number of clusters (rfrnc_qtros) =         52     Root MSE        =   409.6540

                 (Std. Err. adjusted for 52 clusters in countyfip rfrnc_qtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   2.077089    -2.51   0.015    -9.386556   -1.046696
       _cons |    166.374   15.62085    10.65   0.000     135.0139    197.7342
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip rfrnc_
> qtros)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,     51) =       0.95
Statistics robust to heteroskedasticity           Prob > F        =     0.3338
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0005
Number of clusters (rfrnc_qtros) =         52     Root MSE        = 18020.2635

                 (Std. Err. adjusted for 52 clusters in countyfip rfrnc_qtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   81.56063    -0.98   0.334    -243.3277    84.15191
       _cons |  -1395.351   601.9079    -2.32   0.024    -2603.731   -186.9702
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip rfrnc_qtros)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,     51) =       2.12
Statistics robust to heteroskedasticity           Prob > F        =     0.1515
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
Number of clusters (countyfip) =      3,106       Within R-sq.    =     0.0003
Number of clusters (rfrnc_qtros) =         52     Root MSE        = 11151.6317

                 (Std. Err. adjusted for 52 clusters in countyfip rfrnc_qtros)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033   26.58039    -1.46   0.152     -92.0627    14.66203
       _cons |  -1941.697   213.4418    -9.10   0.000    -2370.199   -1513.194
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +------------------------------------------------------------+
     |             var       pm25   earnings   employ~t       lfp |
     |------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217**      -79.6     -38.7 |
  2. | hms_deep_stderr    (0.011)    (2.077)     (81.6)    (26.6) |
  3. |                                                            |
  4. |    observations     75,207    160,346    160,346   161,498 |
     +------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f3.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel f, fourth row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,   1685) =     190.83
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
                                                  Within R-sq.    =     0.0369
Number of clusters (countyfip) =      1,686       Root MSE        =     1.3481

                          (Std. Err. adjusted for 1,686 clusters in countyfip)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0040277    13.81   0.000     .0477402      .06354
       _cons |   9.190593   .0195497   470.12   0.000     9.152249    9.228938
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster countyfip)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   3105) =     123.12
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
                                                  Within R-sq.    =     0.0038
Number of clusters (countyfip) =      3,106       Root MSE        =   409.6540

                          (Std. Err. adjusted for 3,106 clusters in countyfip)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   .4701391   -11.10   0.000    -6.138441   -4.294811
       _cons |    166.374   2.377553    69.98   0.000     161.7123    171.0358
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster countyfip)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,   3105) =      63.44
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
                                                  Within R-sq.    =     0.0005
Number of clusters (countyfip) =      3,106       Root MSE        = 18020.2635

                          (Std. Err. adjusted for 3,106 clusters in countyfip)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   9.991978    -7.97   0.000    -99.17943   -59.99632
       _cons |  -1395.351   50.49432   -27.63   0.000    -1494.356   -1296.345
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,   3105) =      46.47
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
                                                  Within R-sq.    =     0.0003
Number of clusters (countyfip) =      3,106       Root MSE        = 11151.6317

                          (Std. Err. adjusted for 3,106 clusters in countyfip)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033    5.67721    -6.82   0.000     -49.8318   -27.56887
       _cons |  -1941.697   28.67041   -67.72   0.000    -1997.911   -1885.482
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637           0         637     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217***   -79.6***   -38.7*** |
  2. | hms_deep_stderr    (0.004)     (0.470)     (10.0)      (5.7) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f4.dta saved
.                 restore 
.                 
.         }

.         
.         ** panel f, fifth row
.         if 1 {
.                 
.                 * pm2.5
.                 reghdfe pm25 hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster statefip)
(dropped 189 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     75,207
Absorbing 2 HDFE groups                           F(   1,     48) =      18.14
Statistics robust to heteroskedasticity           Prob > F        =     0.0001
                                                  R-squared       =     0.7864
                                                  Adj R-squared   =     0.7634
                                                  Within R-sq.    =     0.0369
Number of clusters (statefip) =         49        Root MSE        =     1.3481

                              (Std. Err. adjusted for 49 clusters in statefip)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |   .0556401   .0130644     4.26   0.000     .0293723    .0819079
       _cons |   9.190593   .0634116   144.94   0.000     9.063096    9.318091
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6677        6677           0    *|
        fe_styr |       637         637           0    *|
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", table(pm25, `tbl_settings_log')
>  replace
(note: file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta not found)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta saved
.                 
.                 * income 
.                 reghdfe d_pc_qwi_payroll hms_deep [aw=seer_pop], a(fe_countyqtroy fe_styr) vce(cluster statefip)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,     48) =      30.23
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.2431
                                                  Adj R-squared   =     0.1760
                                                  Within R-sq.    =     0.0038
Number of clusters (statefip) =         49        Root MSE        =   409.6540

                              (Std. Err. adjusted for 49 clusters in statefip)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -5.216626   .9487319    -5.50   0.000    -7.124179   -3.309073
       _cons |    166.374   4.797857    34.68   0.000     156.7273    176.0208
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631         631           0    *|
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", table(earnings, `tbl_settings_l
> og') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta saved
.                 
.                 * emp 
.                 reghdfe d_pmil_qwi_emptotal hms_deep [aw=seer_pop16plus], a(fe_countyqtroy fe_styr) vce(cluster statefip)
(MWFE estimator converged in 5 iterations)

HDFE Linear regression                            Number of obs   =    160,346
Absorbing 2 HDFE groups                           F(   1,     48) =      11.95
Statistics robust to heteroskedasticity           Prob > F        =     0.0012
                                                  R-squared       =     0.5293
                                                  Adj R-squared   =     0.4875
                                                  Within R-sq.    =     0.0005
Number of clusters (statefip) =         49        Root MSE        = 18020.2635

                              (Std. Err. adjusted for 49 clusters in statefip)
------------------------------------------------------------------------------
             |               Robust
d_pmil_qwi~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -79.58788   23.02065    -3.46   0.001     -125.874   -33.30176
       _cons |  -1395.351   116.3345   -11.99   0.000    -1629.257   -1161.444
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       631         631           0    *|
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", table(employment, `tbl_settings
> _pmil') append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta saved
.                 
.                 * lfp
.                 reghdfe d_pmil_lau_lfp hms_deep [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster statefip)
(MWFE estimator converged in 4 iterations)

HDFE Linear regression                            Number of obs   =    161,498
Absorbing 2 HDFE groups                           F(   1,     48) =      16.75
Statistics robust to heteroskedasticity           Prob > F        =     0.0002
                                                  R-squared       =     0.2403
                                                  Adj R-squared   =     0.1735
                                                  Within R-sq.    =     0.0003
Number of clusters (statefip) =         49        Root MSE        = 11151.6317

                              (Std. Err. adjusted for 49 clusters in statefip)
------------------------------------------------------------------------------
             |               Robust
d_pmil_lau~p |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |  -38.70033   9.454931    -4.09   0.000    -57.71075   -19.68992
       _cons |  -1941.697   47.74823   -40.67   0.000    -2037.701   -1845.692
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |     12424       12424           0    *|
        fe_styr |       637         637           0    *|
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
.                 regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", table(lfp, `tbl_settings_pmil')
>  append
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta saved
.                 
.                 preserve 
.                         use "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", clear
.                         drop if var == "r2"
(1 observation deleted)
.                         drop if var == "_id"
(0 observations deleted)
.                         drop if strpos(var, "_cons_") > 0
(2 observations deleted)
.                         replace var="observations" if var=="N"
(1 real change made)
.                         ingap 3
.                         list , sep(0)

     +--------------------------------------------------------------+
     |             var       pm25    earnings   employ~t        lfp |
     |--------------------------------------------------------------|
  1. |   hms_deep_coef   0.056***   -5.217***   -79.6***   -38.7*** |
  2. | hms_deep_stderr    (0.013)     (0.949)     (23.0)      (9.5) |
  3. |                                                              |
  4. |    observations     75,207     160,346    160,346    161,498 |
     +--------------------------------------------------------------+
.                         saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table1_f5.dta saved
.                 restore 
.                 
.         }

.         
.         
.         
. log close
      name:  <unnamed>
       log:  C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_table1.log
  log type:  text
 closed on:  23 Jun 2022, 23:57:31
----------------------------------------------------------------------------------------------------------------------------------
